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We investigate the symmetric Periodic Anderson Model (PAM) on a three-dimensional cubic lattice 
with nearest-neighbor hopping and hybridization matrix elements. Using Gutzwiller's variational 
method and the Hubbard-III approximation (which corresponds to the exact solution of an appro- 
priate Falicov-Kimball model in infinite dimensions) we demonstrate the existence of a quantum 
critical point at zero temperature. Below a critical value V c of the hybridization (or above a critical 
interaction U c ) the system is an insulator in Gutzwiller's and a semi-metal in Hubbard's approach, 
whereas above V c (below U c ) it behaves like a metal in both approximations. These predictions are 
compared with the density of states of the d- and /-bands calculated from Quantum Monte Carlo and 
NRG calculations. Our conclusion is that the half-filled symmetric PAM contains a metal-semimetal 
transition, not a metal-insulator transition as has been suggested previously. 
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I. INTRODUCTION 

Heavy fermion systems, whose properties are deter- 
mined by nearly localized, strongly correlated /-electrons 
hybridizing with the (i-electrons of the conduction band, 
have been of considerable interest in recent years. □ The 
appropriate theoretical description of heavy fermion sys- 
tems is believed to be the Periodic Anderson Model 
(PAM). Among other properties of heavy fermions, this 
model explains the Kondo effect, i.e., the quenching of 
the magnetic moments of the correlated electrons by the 
conduction electronsH. Recently, using Quantum Monte 
Carlo (QMC) techniques, Huscroft et al.Ha studied a Pe- 
riodic Anderson model with the dispersion of the hy- 
bridization proportional to that of the conduction elec- 
trons. Their study demonstrates that, as the temper- 
ature is lowered, the spin singlet correlation function 
of the conduction electrons develops a sharp structure 
near a critical value of the hybridization. This indicates 
a very rapid cross-over between a Kondo regime and a 
regime where the correlated electrons have unquenched 
moments. The nature of this crossover was also addressed 
by Held and Bullaja who showed that the Periodic An- 
derson model under certain assumptions contains a tran- 
sition equivalent to the Mott-Huhbard metal-insulator 
transition in the Hubbard modclB Since it is clear al- 
ready from Ref. || that one of the assumptions (the strict 
separation of high- and low-energy scales) is at best only 
approximately fulfilled, further analytical and numerical 
studies of this transition in the PAM are clearly called 
for. 

In this paper we study the nature of this 
Mott-Hubbard-like transition analytically, using both 
Gutzwiller's variational methodQ'El and the Hubbard-III 
approximation^. We demonstrate that at half-filling 
there is a quantum critical point as a function of the 
hybridization strength which separates a Kondo regime 
from a phase, in which the d- and /-bands are weakly 



coupled (in the Hubbard-III approximation) or even com- 
pletely decoupled (in Gutzwiller's approach). The disap- 
pearance of the Kondo peak beyond the quantum critical 
point is intimately connected to our choice of the model, 
with a hybridization strength that vanishes at the Fermi 
surface of the conduction electrons. We also present new 
results from 3-dimensional QMC calculations in support 
of the existence of a quantum critical point. 

The grand canonical Hamiltonian describing the two- 
band periodic Anderson model for hybridized d- and /- 
electrons is 

ktr ktr 
i 

icr icr 

Here d Wt7 (f ka ) and d ka (f ka .) are the fermionic operators 
which create and destroy the d-(/-)band electrons of mo- 
mentum k and spin <r, and n i , = f^ a f icr is the number 
operator for the /-electrons of spin a at site i. Further- 
more, Vk is the momentum dependent hybridization term 
between /- and d-electrons. Following Refs. we con- 
sider the dispersion of the d-band and the mixing term 
Vk to be that of nearest neighbor hopping on a three- 
dimensional simple cubic lattice (with unit lattice con- 
stant) whereas the /-band is taken to be dispersionless: 

ek = — 2t [cos k x + cos k y + cos k z ) , (2) 
Vk = — 2P [cos k x + cos k y + cos k z ] , (3) 
e/ = 0. 

Here t and V are the hopping matrix elements between 
the d-d- and /-rf-bands respectively. In this paper we 
study the symmetric PAM in which the chemical poten- 
tial is ji = and (n/) = (n^) = 1. 
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As pointed out also in Refs. HQ there are good rea- 
sons for replacing the usual momentum independent hy- 
bridization, Vw = V, by the momentum dependent hy- 
bridization (g). It follows from elementary symmetry 
arguments that the /- and d- orbitals are essentially or- 
thogonal on the same site. The orthogonality of /- and d- 
orbitals on the same site implies that the hybridization is 
predominantly built up from nearest and further neigh- 
bor contributions. Our choice (||), which assumes only 
nearest neighbor contributions to the hybridization, re- 
flects this fact in the simplest possible manner. As shown 
below, this momentum dependence of the hybridization 
has important consequences, in particular for physical 
properties beyond the quantum critical point. 

This paper is organized as follows. First, in section 
we introduce our two main methods for investigating 
the PAM, namely Gutzwiller's variational method and 
the Hubbard-III approximation. Our variational results 



for the symmetric PAM are presented in section III 



main finding is that the Gutzwiller approach predicts a 
Brinkman-Rice type metal-insulator transition. Next, in 



section IV, we study the PAM in the Hubbard-III ap- 



proximation, which is equivalent to the exact solution 
of a Falicov-Kimball model in infinite dimensions. The 
Hubbard-III solution displays rich behavior as a func- 
tion of the on-site interaction U, including a resonance 
at the Fermi level for weak coupling, a metal-semimetal 
transition at an intermediate-coupling quantum critical 
point, and weakly coupled d- and /-bands at strong cou- 
pling. We then compare the results from the Gutzwiller 
and Hubbard-III approaches to QMC-simulations of the 
PAM on a three-dimensional lattice and also to infinite- 
dimensional QMC-results and to calculations based on 
the Numerical Renormalization Gr oup (NRG, see sec- 
tion |y|) . Finally, in sections VI and |\TJ|, respectively, we 



discuss and summarize our results. 



II. METHODS AND MODELS 

Traditionally, in particular in the context of the Hub- 
bard modeljj there are two famous approaches for investi- 
gating metal-insulator transitions, namely the variational 
approach pioneered by Gutzwilleru and Brinkman and 
Rice,B and the Green function decoupling scheme devel- 
oped by HubbardJj Both of these approaches are clearly 
approximate in nature. Gutzwiller's method predicts the 
formation of an ever narrower quasiparticle peak, accom- 
panied by a divergence of the effective mass, as the on- 
site interaction U approaches a critical value JJ^ utz from 
below. For U > Uf? utz , the Gutzwiller method leads to 
unphysical results, such as the suppression of all hop- 
ping processes and all double occupancies. This method, 
therefore, is more realistic at weak than at strong cou- 
pling. Hubbard's approximation, on the other hand, is 
generally considered to be more realistic at strong cou- 
pling. At weak coupling it predicts the steady decrease of 



the number of charge carriers at the Fermi level; however, 
this mechanism is implemented in such a way that Fermi 
liquid properties are violated. At strong coupling Hub- 
bard's method predicts band splitting, i.e., the formation 
of a lower and an upper Hubbard band. Combination 
of both methods (the Gutzwiller method at weak and 
Hubbard's approximation at strong coupling) has yielded 
valuable information on the metal-insulator transition in 
the Hubbard model. Here we combine both methods in 
order to shed light on the nature of the metal-insulator 
transition in the PAM. 

To study the ground state properties of the half-filled 
PAM, Eq. (|IJ), within the Gutzwiller approach, ra^fo!-. 
low the variational procedure of Rice and UedaJJiijilII 
These authors considered the PAM with on-site (rather 
than nearest-neighbor) hybridization. The central aspect 
of the Rice-Ueda approach is the suppression of doubly 
occupied /-states. The Gutzwiller-correlated wave func- 
tion, \iPg)i is f° r the case °f the PAM defined as 



\1> a ) = P\ih) , 



(4) 



where \ tpo) is the wave- function for the uncorrelated (U = 
0) ground state at half-filling, and P is the Gutzwiller 
correlator, defined as 



P = g d =l[[l-(l-g)Di 



(5) 



The operator D = Yli ri if1 n ifl m Eq. (|B|) is the double 
occupancy operator for the /-electrons and g is a vari- 
ational parameter. For g = 0, the operator P projects 
all the states onto the subspace without doubly occupied 
/-sites, whereas 5 = 1 corresponds to the uncorrelated 
state. In general g has to be determined by minimizing 
the total energy of the system in the correlated state. 

For the symmetric PAM, which is the case of inter- 
est in this paper, the Gutzwiller correlator P is treated 
by renormalizing all hopping processes by a Gutzwiller 
factor q(d), where d = D/Af is the fraction of dou- 
bly occupied /-sites and Af is the total number of lat- 
tice sites. This approximation, which is alternatively re- 
ferred to as the "Gutzwiller approximation" or as "semi- 
classical counting", becomes exacmA (at least within 
the Gutzwiller variational approach, not for the PAM- 
Hamiltonian itself) in the limit of high spatial dimensions 
(d = oo). 

The central element in Hubbard's Green function de- 
coupling scheme is the so-called "alloy analogy" , in which 
it is assumed that the down-spins hop while the up-spins 
are immobile, and vice versa.Ej The Hubbard-III approx- 
imation, like the Gutzwiller approach, can be understood 
as the exact solution of a simplified problem in high spa- 
tial dimensions. In the context of the PAM, the mobile 
nature of one spin species (say the /^ -spins), interacting 
with an "alloy" of immobile electrons of opposite spin 
(here the /j -spins), can be described by the following 
Hamiltonian: 
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H = -tY,dtd ja -vj2(4 l fu+Mi) 



u 



i/- 2) [ n ifi 2 



E e / n i/--^E("i/- +n i^) 



In this simplified model, the /^-electrons form an alloy 
of immobile spins, since their hybridization with the ri- 
band vanishes. As a consequence, the dj-electrons are 
completely decoupled from the rest of the system and can 
be integrated out. The creation (annihilation) operators 
for the remaining (^-electrons will simply be denoted by 
d\ (rfj ) below. Since we consider the symmetric PAM 
(e/ =0) at half- filling (fi — 0), the last two terms in the 
Hamiltonian drop out, and we are left with 



,d,' 



H = -tj24d i -vj2(4f il ,- Jirj , 

(ij) (ij) 

'E( n i/T~2-) ( n ifl-h) ■ ( 6 ) 
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This model is very similar to the spinless two-band 
Falicov-Kimball model, for which the Hubbard-III ap- 
proximation (in the framework of the Hubbard model) 
becomes exact in the limit d — > 00. Analogously, the 
appropriate Green function decoupling ( "alloy analogy" ) 
for the PAM becomes exact in the simplified model (||) 
in high dimensions. Below, in section IV, we solve the 
model (^) exactly in d = 00. In order to obtain meaning- 
ful results in this limit, the hopping and the hybridization 
have to be scaled as t = i/^/Z and V = V/V~Z, respec- 
tively, where Z denotes the number of nearest neighbors 
of each lattice site. 

Part of the Hubbard-HI approximation is the choice 
of the lattice. Hubbardta started from a semi-elliptical 
density of states (DOS), which corresponds to the Bethe 
lattice in high dimensions. Below, we follow Ref. |j] and 
solve (H) on the Bethe lattice. This choice has several 
advantages, e.g., that the bandwidth is finite in d = 00, 
that the DOS near the band edges resembles the DOS on 
a simple cubic lattice in d = 3, and that the properties 
of the DOS can be studied analytically, since one obtains 
a relatively simple closed equation for the local Green 
function. 

Essentially the same Hamiltonian (|J), but now with 
on-site hybridization, was considered also by Consiglio 
and Gusmao.Ej These authors referred to the model as 
the Simplified Periodic Anderson Model (SPAM), a des- 
ignation that we extend also to the case of more gen- 
eral (in particular nearest-neighbor) hybridization. The 
method [Of solution in Ref. [14| was that of Bijandt and 
Mielsch.Ej Here we use an alternative methodtj that is 
much better suited for the calculation of the DOS on the 
Bethe lattice. An early application of the Hubbard-III 
scheme to the periodic Anderson model, in particular a 



calculation of the resistivity as a function of temperature, 
can be found in Ref. [It]. 



III. GUTZWILLER'S VARIATIONAL APPROACH 

In the Gutzwillcr approximation, the kinetic energy 
terms of the PAM are mapped to an effective Hamil- 
tonian with a renormalized hybridization. This result 
was first obtained with the use .af.-semiclassical count- 
ing arguments by Rice and UedaEltll (see also Ref. fl8| ), 
and then put on a solid footing by Gebhardjl3 who 
showed that this so-called "Gutzwiller approximation" 
becomes exact in the limit of infinite spatial dimension- 
ality We recall that Gutzwiller's variational scheme is 
equivalent to the-slave-boson mean-field theory of Kotliar 
and Ruckensteintj at T = 0. At half-filling, the efeptive 
Hamiltonian in the Gutzwiller approach becomes! 1 3uJ 

H cS = 'A***, + E ^ (4.4. + flM , (?) 

kcr kcr 

with the renormalized hybridization 



Here the renormalization factor q(d) takes the forn: 



q = 8(2(1 - 2d) 

and d = D/N is the fraction of doubly occupied /-sites. 
The ground state energy E g is now obtained from the 
expectation value of H e s, 

E g = {MH cS \^)/N + Ud- U/4 . 

The effective Hamiltonian is easily diagonalized. For con- 
venience, we set t = 1, to establish a unit of energy. One 
then finds two bands, with eigenenergies i?±ek, where 
i9± = |(1 ± ^1 + 4qV 2 ). Note that there is no gap be- 
tween the two bands, and that both bands contribute to 
the ground state energy. This is in .marked contrast to 
the results for on-site hybridization£3'U The result for 
the ground state energy is: 



E g = — |e |y 1 + 32d(l - 2d)V 2 + Ud-U/A, (8) 

where eo = X){ £k <o o-} e ko-/A/' is the ground state energy 
per lattice site of a completely decoupled d-band (U — 
V — 0) . Minimizing the ground state energy with respect 
to d yields 



(9) 



Here we introduced the critical value U c of the interaction 
or, equivalently, a critical value V c of the hybridization, 
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FIG. 1. The fraction of doubly occupied sites is plotted as 
a function of 1/V, where V is the hopping between the /- 
and the d-bands. For the critical V c , given in Eq. @, the 
number of doubly occupied sites vanishes and all the sites 
become localized. The parameters are chosen to be t — 1 and 
U — 6, and the average energy of the system is calculated by 
summing over all the filled bands: |eo| = 60/4 3 , where the 
lattice size is taken to be 4 x 4 x 4. With the above values of 
the parameters the critical value of V is approximately 0.63. 



for which the expectation value (^) of the double occu- 
pancy vanishes: 



U c = 16V 2 \e \ 



V c 



u 
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(10) 



only applies for U < U c (or V > V- 



Since d > 0, Eq. j omy ii ppnes u»r < •_. < (or i • 
For U > U c , d = 0. Substituting Eq. (§) into Eq. 
gives the ground-state energy of the system: 



E g = -\e Q W[l + (2V) 2 ] [l + (2UV/U c y 



(11) 



We note that the variational ground-state energy is inde- 
pendent of the hybridization, E g — — ^ — |eo|, for U > U c . 



In Fig. |l| we plot the fraction of doubly occupied /-sites 
d as a function of the inverse d-/-band hopping parameter 
V. The fraction d decreases monotonically with V and 
finally vanishes at a critical value V c . This implies that, 
at this critical value of V, every /-site is singly occupied, 
i.e., that the /-electrons are localized. 

The critical hybridization strength V c (or, equivalcntly, 
the critical interaction U c ) defines a quantum critical 
point which separates two distinct regimes in the model. 
For V > V c (or U < U c ), the /-electrons hybridize with 
the d-electrons and their moments are screened. This is 
the Kondo regime. On the other hand, for V < V c (or 
U > U c ), the two bands are decoupled at the Fermi sur- 
face, and the /-electron moments are unquenched. Thus 
the /-electrons show a "metal-insulator transition" from 
extended (U < U c ) to localized (U > U c ). In contrast, 



the d-electrons at the Fermi level are always extended in 
both regimes. ■— ■ 

Following Rice and Uedajl^l we introduce the binding 
energy in the Kondo regime as the singular part of the 
variational ground state energy, 

E b (U) = -U/4-\e \-E g (U) , 

so that E b (U) = for all U > U c . With this definition of 
the binding energy, it is easy to show that E b is positive 
for all U < U c , and that for U \U C (i.e., if the transition 
is approached from below): 



E b 



2y 2 M 
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Near the transition, where the two energy scales E b and 
| eo j are well-separated, it seems plausible that E b can 
be identified with the Kondo temperature, E b ~ /crIk- 
It is of interest Jxi compare our result for Ik to that of 
Rice and Ueda,Q who found T K = 2\e \e- u ^ 8v2 . While 
these results cannot be compared in detail (since Ref. 

assumes a one-dimensional, linear dispersion for the 
d-electrons) , it is nevertheless clear that the Kondo tem- 
perature in our model is strongly suppressed relative to 
that of Ref. [H] for U < U c . 

The density of states of the PAM in the Gutzwiller 
approach is of interest, too. If we denote the DOS of a 
decoupled d-band (U = V = 0) by v{E), then the DOS 
for the interacting d-electrons in the Gutzwiller approxi- 
mation is simply given by 
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Thus one finds that the DOS of the ci-electrons near 
the transition (U < U c ) is strongly enhanced at the 
Fermi level compared to the situation for U > U c , where 
Vd(E) = v{E). This result represents an interesting 
counter-example to the exhaustion scenario of Nozieres.tj 
At low temperatures (T < Tk) only the electrons within 
7k of the Fermi surface can effectively participate in 
screening the local moments. In a concentrated system 
such as ours, there are more moments to screen than 
conduction-band states available for screening. In the 
metallic regime of the conventional PAM, this should lead 
to a depletionj-Qf— the density of screening states at the 
Fermi surface.rHr 2 ! However, within the Gutzwiller ap- 
proximation this "exhaustion physics" is clearly absent, 
as we see an enhancement of the (i-band DOS near the 
Fermi surface. 

Similarly the /-band DOS takes the form 



ME) 



(tf + /|0-IMW-l) + (W-\/^ + HE/^ 
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Thus, in the Gutzwiller approximation, the quasiparticlc 
peak has width oc q oc Jcx (1-U/U C ) oc (|e |T K ) 1/2 , 
which is much larger than the Kondo scale Tk. Hence, 
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interestingly, the Kondo temperature Tk near the tran- 
sition is not determined by the width of the peak in the 
/-DOS, but rather by the (much smaller) binding energy 
Ei,. The physical explanation for this is that screening 
becomes increasingly less efficient as one approaches the 
transition, due to the renormalization of the hybridiza- 
tion rates, Vk = y^^k — > for U — ► U c . 



IV. THE HUBBARD-III APPROACH 

The Hubbard-III approach (alloy analogy) is equiva- 
lent to the exact solution in d — oo of the Simplified Pe- 
riodic Anderson Model (SPAM), Eq. (|) with t = t/VZ 
and V = V/yfZ. For the calculation of the DOS on the 
Bethe lattice, it is most convenient to first map the SPAM 
to an effective non-interacting Hamiltonian, following the 
lines of Rcf. 16, aad then use the renormalized pertur- 
bation expansionE3 to calculate the DOS. This method 
also allows one to conclude immediately that the DOS at 
half- filling is temperature independent, as a consequence 
of the non-interacting nature of the effective Hamilto- 
nian. Along the lines of Ref. |l^ we find that the Fourier 
transform of the local matrix Green function, 



Gii(r) 



(rd ; (r)4(0)> (Td,(r)/t(0)) 
(T/ U (r)4(0)) (T/ u (r)/t(0)) 



satisfies a cubic matrix equation, 

G{z) = ± {[/_ - ©GOf 1 + [/+ - OGOr 1 } , ( 
with 
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In spite of the formally simple structure of Eq. [lj, the 
detailed analysis of the DOS of the d- and /-elecfapns is 
rather involved and will be published elsewhere.Eil Here 
we focus on the physical content of the Hubbard-III re- 
sults and compare them to results from the Gutzwiller 
approach and various numerical techniques. For conve- 
nience we put t — 1 to fix the unit of energy; note that 
this convention differs from that of the previous section. 



A. Density of states at strong interaction 

First we consider the results from the Hubbard-III ap- 
proximation in the limit of large interaction, U — > oo. In 
this limit almost all the spectral weight of the d-electrons 
is contained in a semi-elliptic band near E = 0: 



Z7T 



(13) 



Similarly, nearly all spectral weight of the /^-electrons 
is contained in two high and narrow peaks near E = 



±i£7, whose width is of order V 2 /U while their height is 
proportional to U/V 2 . With a redefinition of the energy 
variable as A = U(\E\ - ^U)/V 2 , one finds 



u 



4ttV 2 \ 



V6A - 1 - A 2 



provided that the argument of the square root is pos- 
itive (3-\/8<A<3 + \/8); otherwise v fi (E) van- 
ishes. The presence of an upper and a lower Hubbard 
band in the fi -spectrum at first sight suggests the occur- 
rence of a metal-insulator transition at some finite value 
of U. However, closer inspection shows that there is some 
small additional spectral weight for the ci-electrons near 
E = ±±U, 
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and small additional spectral weight for the /^-electrons 
near E — 0: 



2V 2 
ttU 2 



E 2 V^~E 2 



This is a first clear indication that the /j -electrons in 
the SPAM do not undergo a metal-insulator transition. 
Instead one finds a metal- semimetal transition: at arbi- 
trarily large U, the d- and /^-bands remain weakly hy- 
bridized. 



B. Density of states at the Fermi level 

Exactly at the Fermi level (E — 0) an explicit non- 
perturbative solution can be obtained for all interaction 
strengths U. A detailed analysis of the cubic matrix 
equation (|l^) shows that there is only one physically ac- 
ceptable solution, 



G(0) = 



fli(0) 52(0) 
52(0) 53 (0) 



where the matrix elements of G(0) are given by 
5 i (0) = -V2M+ 



and 



52(0) 
53(0) 



R± 



V2V 



V2V- 



(R+-R-) 
(R+ - R-) 



U 2 



ill 



u 2 



These results hold only for U < C/ AA = 2V 2 , where 
the superscript "AA" stands for "alloy analogy". For 
all U > L/ AA the solution is simply given by g\ (0) = — i 
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and 52(0) = g 3 (0) = 0. Note that the Hubbard-III and 
Gutzwiller approaches predict the same dependence of 
U c on the hybridization: both U AA and U^ utz are sim- 
ply proportional to V 2 . 

Physically these results mean, that at weak coupling 
the DOS at the Fermi level decreases fairly slowly as the 
interaction U is switched on, both for the d- and the 
fl -electrons: 



v sm ~ ^ [}- w*) 

while the critical values are approached quite rapidly: 



Mo) - ~ [i + yu AA -u/v 



(u t u AA ) . 



For all U > U AA the /|-DOS at the Fermi level van- 
ishes exactly, while the <i-DOS is pinned at the value 
^d(O) = 1/tt. This clearly demonstrates that U AA marks 
a quantum critical point, although the nature of the tran- 
sition cannot be deduced from an investigation of the 
DOS only at the Fermi level. 

We add that, at weak coupling, the DOS for both the 
d- and the /-electrons displays an interesting resonance 
at the Fermi level. This resonance is very narrow, of 
0(U 2 ), and its height remains of 0(1) for U ->■ 0. The 
amplitude of the resonance for the d-electrons is positive, 
so that the (i-band DOS at U = + is larger than for 
[7 = 0. In contrast, the amplitude of the resonance in 
the /-band is negative, i.e., the /-band DOS at U = + 
is smaller than for [7 = 0. 



C. Density of states at and beyond the critical point 

We now consider the shape of the DOS as a function of 
energy, first for U = U AA and then for U > U AA . We fo- 
cus on the energy interval near the Fermi level, since this 
interval determines the nature of the quantum critical 
point and virtually all physical properties of interest. 

First we present the result for the DOS exactly at the 
transition, i.e., for U — U AA . In this case one finds 
that the DOS has a sharp (cubic-root) singularity as a 
function of energy for \E\ — > 0: 



^(£)~i + #(l + ^) 1/3 |£| 1/3 



1 ^ 1/3 |R|V3 



(|E|-0) 



Cubic-root critical behavior of the DOS near the Fermi 
level is well-known from the Hubbard-III solution for the 
Hubbard model. Here we find that similar critical be- 
havior occurs also for the SPAM with nearest-neighbor 
hybridization. 



For energies near the Fermi level, one generally finds 
that the d-DOS for U > U AA is metallic while the /j- 
DOS is semimetallic. In particular, one finds that the 
/j_-DOS contains spectral weight, arbitrarily close to the 
Fermi level, 



w 2 
rtf 2 



E 



(\E\^0,U>U ( 



Note that the coefficient in front of the semimetallic E 2 - 
behavior is valid for all U > U AA and, h ence, simply 
identical to that calculated in section IV A for large U. 



V. COMPARISON TO QMC AND NRG RESULTS 

To obtain a complete picture of our model, we com- 
pare our analytical results to single-particle spectra ob- 
tained from three different numerical techniques: (i) 
QMC simulations of the PAM on a three-dimensional lat- 
tice (where present new results and compare to results 
from recentB'Q calculations), (it) QMC .simulations of the 
PAM on an infinite-dimensional latticea, and (Hi) recent 
results from NRG calculations in infinite dimensions!} 

We start with the comparison to the three-dimensional 
QMC-results. In Fig. ||we plot new results for the <i-band 
DOS obtained using the maximum entropy method to 
analytically continue three-dimensional QMC data calcu- 
lated at T = 0.2. For small V the presence of upper and 
lower Hubbard bands well separated by a gap U can be 
seen in the /-band DOS. With increasing V, the weight 
in the central region is enhanced at the expense of the 
Hubbard side-bands. 

For V larger than the critical hybridization V c , the 
DOS for the <i-band is greatly enhanced compared to the 
non-interacting (i-band DOS (shown in dots). This is 
consistent with the previously mentioned enhancement of 
the d-band DOS in the Gutzwiller approximation in this 
regime. On the other hand, in the regime where V < V c , 
the Green function is not significantly different from the 
non-interacting Green function. This indicates that for 
V < V c , the d-band is weakly coupled to the /-band. 
This result is compatible with both the Gutzwiller result 
(that the d-band is strictly ^-independent for V < V c , 
i.e., v d (E) = v(E)) and the Hubbard-III result, Eq. 
which holds asymptotically both for U — > 00 and V — > 0. 

The /-band DOS was studied in Refs. |J| with the 
use of finite- and infinite-dimensional QMC techniques. 
The infinite-dimensional (DMFT) calculationso were also 
done for a strictly paramagnetic state, without antiferro- 
magnetic fluctuations. For small hybridization, V < V c , 
one finds that the /-band DOS consists primarily of an 
upper and a lower Hubbard band around E = ±5 £7 7 re- 
spectively. Only in the DMFT calculations^ is there a 
sign of small spectral weigth near the Fermi level. For 
increasing V , one finds that the upper and lower Hub- 
bard bands merge. Before the bands merge (V < V c ), 
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they also point out that one of these assumptions (the 
strict separation of high- and low-energy scales) is at 
best only approximately fulfilled for U > U c , especially 
for larger values of the hybridization rate (V > 1). The 
Hubbard-III results of this paper strongly suggest that 
the high- and low-energy scales in the PAM are in fact 
never strictly separated, so that the transition is from 
a metal to a semimetal, not an insulator. Physically it 
seems obvious that for large U, as a result of the finite 
hybridization rate, the /-band must have some spectral 
weight near the Fermi level, where the weight of the ri- 
band is concentrated. Similarly, the ci-band must have 
some spectral weight near the upper and lower Hubbard 
bands of the /-electrons, i.e., near E — ±\U . 
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VI. DISCUSSION 



FIG. 2. The d-band density of states, ua(fj), for the 3d 
PAM via finite-d QMC for various hybridizations, V > V c . 
From Refs. | and | we estimate that V c ~ 0.5. The DOS at the 
Fermi energy is resonantly enhanced over the non-interacting 
value, shown as dots in the figure. For V < V c , the Green 
function is similar to the non-interacting Green function, 
which indicates that the d- and the /-bands are approximately 
decoupled in this regime (not shown). 

additional resonant peaks develop, which are character- 
istic for Kondo singlet formation. After the bands merge 
(V > V c ), the /-band DOS in the finite-dimensional sim- 
ulation shows some depletion near the Fermi level. Since 
this depletion is absent in the DMFT calculations, this 
feature is most likely due to short-range antiferromag- 
netic fluctuations. Comparing these numerical results to 
those from the Hubbard-III and Gutzwiller approaches, 
we observe that the merging of the two bands is described 
by the Hubbard-III approach but not by the Gutzwiller 
approximation. The Hubbard-III approach also explains 
the small additional spectral weigth near the Fermi level 
in the DMFT calculation. The Gutzwiller approxima- 
tion cannot explain cither the upper and lower Hubbard 
bands or small spectral weigth near the Fermi level. Ob- 
viously, neither approach can explain the depletion due to 
short-range antiferromagnetic fluctuations in the finite- 
dimensional calculations. n 

We add that, very recently, Held and BullaO per- 
formed calculations for the paramagnetic ground state 
of the PAM using the Numerical Renormalization Group 
(NRG). These calculations clearly reveal significant ad- 
ditional spectral weight near the Fermi level for U > U c . 
The numerical results in Ref. ^ are the clearest evidence 
to date that the quantum critical point in the PAM corre- 
sponds to a metal-semimetal and not to a metal-insulator 
transition. Interestingly, the authors of Ref. || show an- 
alytically that the PAM under certain assumptions con- 
tains a transition equivalent to the Mott- Hubbard metal- 
insulator transition in the Hubbard model. However, 



The two main questions to be discussed are: What can 
be learned from the Gutzwiller and Hubbard-III approx- 
imations for the PAM? And what features are missing in 
these two analytical approximation schemes? 

As pointed out before, in section H one expects 
Gutzwiller's variational method to be more realistic at 
weak than at strong coupling. This expectation is based 
on experience with the Hubbard model and on the fact 
that the Gutzwiller wave function becomes exact in the 
weak coupling limit. For similar reasons, Hubbard's ap- 
proximation is considered to be more realistic at strong 
coupling. For example, the Hubbard-III approximation 
violates weak-coupling (Fermi-liquid) properties and is, 
therefore, better suited for U > U c , where the DOS of 
the /-electipirs displays essentially no low-energy spec- 
tral weight £a Combination of these two approaches leads 
to the following physical picture. 

We distinguish high and low spatial dimensions. In 
high spatial dimensions, short-range antiferromagnetic 
fluctuations (in particular due to spin-flip processes on 
nearest-neighbor sites) are small, which is a prerequisite 
for the validity of the Gutzwiller and Hubbard-III ap- 
proaches. The Gutzwiller approach then shows that the 
/-band DOS develops a quasiparticle peak near the Fermi 
level if the interaction U is turned on. The quasiparti- 
cle peak becomes ever narrower, until it vanishes at the 
critical interaction U c , where the metal turns into a semi- 
metal. The Gutzwiller approach is unable to describe the 
side bands in the DOS that develop for U < U c . These 
high-energy features are better captured by the Hubbard- 
III approximation which, in turn, is unable to describe 
the Fermi-liquid peak. For U > U c one expects the DOS 
to remain semi-metallic, as is seen in the Hubbard-III 
approach. The Hubbard-III results further suggest that, 
at some finite interaction U > U c , the DOS splits into 
three parts: a lower Hubbard band near E — —hU, a 
semimetallic low-energy part of the spectrum, and an up- 
per Hubbard band near E = +§Z7. This scenario is con- 
sistent with the results from the high-dimensional QMC 
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simulations and agrees very well with those of the NRG 
calculations discussed in section |v|. 

In contrast, there will be significant short-range an- 
tiferromagnetic fluctuations in the paramagnetic phase 
in low spatial dimensions (certainly for d = 1,2 and 
to some extent also for d = 3). In these dimensions 
one, therefore, expects depletion of the /-band DOS 
near the Fermi level. This depletion is..explicitly seen 
in three-dimensional QMC-calculationsBa These deple- 
tion effects cannot be described in terms of the param- 
agnetic Gutzwiller wave function for the PAM, since it 
does not contain the necessary antifcrromagnetic fluctu- 
ations. The effects of finite-dimensionality on the quality 
of the Hubbard-III results for U > U c is probably not 
large. One expects the semimetallic part of the DOS to 
be somewhat depleted in d = 3 as compared to d = oo. 
An interesting question that cannot be answered at this 
stage is, whether the semimetallic pseudogap changes 
its functional form due to short-range antifcrromagnetic 
fluctuations and, in particular, whether these fluctua- 
tions might turn the pseudogap into a hard gap. Given 
these uncertainties concerning the influence of antifcrro- 
magnetic fluctuations in finite dimensions, it seems fair 
to say that our results are at least consistent with the 
three-dimensional QMC data published here and in Refs. 
H and ||. We do not believe that our results can reason- 
ably be applied to one- and two-dimensional systems, 
since physics in these low dimensions differs too much 
from that in high dimensions, where our approximations 
are valid. 

In the Gutzwiller and Hubbard-III methods studied in 
this paper, only states without broken symmetry were 
considered. The true ground state of the half-filled sym- 
metric PA M m pfj course , may well be antiferromagnetically 
ordered. QoEjO In this sense our results are not so much 
relevant for the ground state of the PAM, but rather 
for the paramagnetic phase at slightly elevated temper- 
atures (T > TNeei)- Note that the low-temperature an- 
tifcrromagnetic phase can Jae partially or entirely sup- 
pressed due to frustration^ With this proviso, we be- 
lieve that the main result of our paper, the absence of 
Kondo screening for V < V c (or U > U c ), is robust. 
In particular, we expect that metal-semimetal transition, 
predicted by the Hubbard-III approach, actually occurs 
in the PAM with nearest-neighbor hybridization. 

A comment is in order concerning the critical val- 
ues of the interaction in the Gutzwiller and Hubbard- 
III ("alloy analogy") approaches, [/f utz and U AA . In 
sections |ffl| and [TV] we found that £/ c Gutz = 16V 2 \e \/t 2 
and U AA = 2V 2 /t, where we reinstated factors of t and 
t = i/\/~Z. In order to be able to compare these two re- 
sults, we calculate the critical value [/r utz for the Bethe 
lattice with coordination number Z — > oo. One read- 
ily finds that in this case | eo j = 8i/37r, so that U~ iutz — 
128V 2 / (37rf). Hence there is a significant discrepancy be- 
tween the two [/ c -values: t/ Gutz /t/ AA = 64/3tt ~ 6.79. 
Partly this is due to the neglect of "resonance broad- 



ening corrections" , which are part of the full Hubbard- 
III approximation but aee not taken into account in the 
alloy-analogy approachJla For the Hubbard, model, the 
resonance broadening corrections are knownQ to enhance 
U AA by a factor of \/3, so that one then finds a ra- 
tio f/ c Gutz /(V3C/ AA ) = 32/3\/37r ~ 1.96. Hence the 
Gutzwiller prediction for U c in the Hubbard model is sig- 
nificantly larger than the full Hubbard-III result. In the 
PAM the discrepancy is even larger. Even if resonance 
broadening corrections are taken into account through an 
ad hoc factor of %/3, t/ Gutz is still larger than V3U AA by 
a factor of nearly 4. In order to determine which of the 
three predictions £/ Gutz , U AA and \/3t/ AA is best, we 
compare with results from recent NRG-calculations for 
the ground state of the PAM in high dimensions.El E.g., 
for i = V = 1 one finds that U^ AM ~ 4. In conclusion, 
compared to the NRG- value, U AA is too small by a factor 
of 2 due to the neglect of resonance broadening correc- 
tions, the "corrected" value ^/3U AA is quite close, and 
the Gutzwiller prediction J7 Gutz is too large by a factor 
of ~ 3.395. For the PAM, therefore, the full Hubbard-III 
approximation leads to a much more accurate prediction 
for the critical interaction than the Gutzwiller approach. 

We now comment on the influence of small devi- 
ations from the strict nearest-neighbor hybridization, 
V£ s = ~2VJ2 ( LiCOs(k e ), which we considered through- 
out in this paper. As an example, we consider combi- 
nations of this "extended s-wave" form and a local ( "s- 
wave") hybridization, V£ — —V s , or a "p-wave" form, 
V£ = — 2V P X^=i sin(fcf). In the Gutzwiller approxima- 
tion, one readily finds that there is now a gap at the 
Fermi surface and that the renormalization factor q is 
now finite for all U > 0. For instance, one finds for a 
mixture of an extended s-wave and a s-wave: 



qoc^exp 



U - U C (V) 



Wv(0)V 2 



(U > U c ) 



(14) 



For a mixture of an extended s-wave and a p-wave, the 
factor V 2 in the exponent is replaced by the average of 
(V^T) over the non-interacting Fermi surface of the d- 
electrons. Strictly speaking, Eq. ( |l4] ) implies that the 
quantum critical point, which occurs for V s = V p = 0, 
is unstable with respect to small s- or p-wave pertur- 
bations. However, Eq. ( p^ ) also shows that the QCP at 
U C (V) = 16F 2 |eo|/£ 2 is replaced by a sharp crossover to a 
state with very heavy quasiparticles. Experimentally this 
could hardly be distinguished from the metal-insulator 
(or metal-semimetal) transition found for V s = V p = 0. 

Finally we comment on the difference between the local 
and nearest-neighbor hybridization rates from a renor- 
malization group point of view. The Kondo effect is due 
to the antiferromagnetic coupling of d-electrons near the 
Fermi surface to the /-electron moments. As the system 
renormalizes, the coupling scales to the strong coupling 
(J = oo) fixed point. Adjacent to this fixed point, only 
energy scales and momenta near the Fermi surface are 
important. For the standard PAM the hybridization is 
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constant, and one finds the usual Kondo effect. In con- 
trast, the most salient feature of the model we study is 
that the f-d hybridization vanishes at the Fermi surface 
of the half filled non-interacting model, as can be seen 
from Eqs. (^) and (^). The corresponding exchange cou- 
pling also vanishes on the Fermi surface, and therefore 
cannot be rescaled to the strong coupling limit. Thus 
the absence of the Kondo effect in our model at small 
values of V/t is fully consistent with the usual renormal- 
ization group arguments. 



VII. SUMMARY 

We have studied the half-filled symmetric Periodic An- 
derson model with an f-d hybridization proportional 
to the (i-band dispersion, using Gutzwiller's variational 
wave function and Hubbard's alloy-analogy approach. 
Both methods demonstrate the occurrence of a quan- 
tum phase transition for a critical value V c of the f-d 
hybridization rate or, equivalently, for a critical interac- 
tion strength U c . For V > V c (or U < U c ) the system 
is a metal. For V < V c (or U > U c ) the system be- 
haves as an insulator in Gutzwiller's approximation and 
as a semi-metal in the Hubbard-III approach. Based on 
the results of these two approximate methods, we predict 
the occurrence of a similar quantum phase transition in 
the exact solution of the PAM with nearest-neighbor hy- 
bridization. A number of the properties of this transition, 
such as a strong suppression of the Kondo temperature 
on the metallic side, were also predicted. Most impor- 
tantly, we conclude from our results that the nature of 
the transition in the exact solution of the PAM .will be 
that of a metal- semimetal, not a metal-insulatora tran- 
sition. These analytical results are compared with the 
DOS obtained from QMC-simulations and NRG-results 
for the d- and /-bands. We find good agreement between 
our scenario and QMC- or NRG-results in high spatial 
dimensions (d = oo), while our results are at least con- 
sistent with the QMC-calculations published here and in 
Refs. H and || for the three-dimensional PAM. 
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